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ABSTRACT 


We consider a class of grand unified theories in which cosmologi cal ly 
significant axion and neutrino energy densities arise naturally. To obtain 
large scale structure we consider (1) an inflationary scenario, (2) inflation 
followed by string production, and (3) a non-infl ationary scenario with 
density fluctuations caused solely by strings. We show that inflation may be 
compatible with the recent observational indications that !2 < 1 on the scale 
of superclusters, particularly if strings are present. 
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Axions with mass on the order of 10" 3 - 10‘ 4 eV, have been suggested as 
candidates for the dark matter in galactic halos 1 * 2 . It has also been shown 
that axions with a cosmological ly significant energy density provide an 
important component in the mechanism for generating structure in the universe 
on scales up to 10 15 Mq 3 » 4 . In this picture, axions, being gravitationally 
unstable on all scales, will cluster first, providing the seed potential wells 
for galaxy formation so that the galaxy distribution on scales up to ~ 10 15 Mq 
clusters would naturally follow the axion mass distribution. Observational 
support for such a relationship is discussed by Blumenthal et al . 5 . They 
point out that the ratio of dark to luminous mass is roughly constant up to 
the scale of rich galaxy clusters. 

J* 

An S0(10) GUT framework which leads to the production of cosmological ly 
significant axions has been given 6 . In this letter, we first argue that 
within this class of models (and suitable extensions thereof such as Eg), a 
cosmological ly significant neutrino mass is obtained naturally. We then 
proceed to discuss some cosmological implications of this result for the 
formation of structure in the universe within the context of three different 
scenarios, (1) an inflationary scenario, (2) an inflationary scenario followed 
by string production, and (3) a non-inflationary scenario with density 
fluctuations produced solely by strings. 

As an example of a grand unified theory which gives n - 8 , consider 
the following S0(10) model 6 (the global U(l) Peccei -Quinn symmetry 7 is not 
explicitly exhibited): 

S0(10) ► SU ( 3 ) x SU( 2) - x SU(2 ) d x U(1) r , 

M x ~10 15 GeV L K B ' L 

► SU(3) x SU(2). x U( 1) ► SU(3) x U(l) 

f - 10 12 GeV L M w ~100 GeV em 


( 1 ) 
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Both the global U ( 1 ) symmetry and the local B-L symmetry are broken at a 
scale of order 10 12 GeV. (Note that the value of the intermediate scale is 
not put in by hand, but is determined from the renormalization group equations 
of the gauge couplings). From the results of Reference (1), it follows 
that ii * 0.1-1. 

d 

Let us now consider neutrino masses in this model. The breaking of B-L 
at scale f a , caused by a ^ -plet of Higgs fields, induces a Majorana mass 
term for the right-handed neutrino v Rl of order h^f a , where h-j denotes the 
Yukawa coupling of the 1 th generation. The breaking of SU(2) x U(l) to U(l) em 
is achieved by a Higgs 10 plet and gives rise to Dirac mass terms m^) = m 

i u i 

(where u-, denotes u,c,t,...) linking the left and right-handed neutrinos. 
Moreover, it can be shown that an effective Majorana mass term for the left- 
handed neutrino v of order c^ - h 1 (X^/X^) << ^io > ^^a a ^ so ^ nclucec|8 * Here 
Xj denotes the quart! c higgs coupling between the ^£6 and the ^Q, X ^ is the 
quartic self-coupling of and < 4 > jg > is the vacuum expectation value of 

the LQ,. With f * 10^ GeV, X,/X„ of order unity, and h. - O(g^) (where g 
denotes the S0(10) gauge coupling), c 1 is in the electron volt range. 

Di agonal ization of the neutrino mass matrix (neglecting, for simplicity, 
mixings between generations) yields the eigenvalues 


^heavy - ^l^a* 

("’flight c i "u/KAeavy* 


( 2 ) 


It follows from eq. (2) that electron volt neutrino masses arise 
naturally in the class of models under discussion. Indeed, due to the 
presence of the c 1 term in the mass matrix, the light neutrino of each 
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generation can have a mass in the electron volt range. Thus, neutrinos can 
contribute significantly to the dark matter in the universe. 

We now discuss the implications of significant axion and neutrino energy 
densities for the evolution of structure in the universe. Two mechanisms for 
producing density fluctuations in the early universe have been extensively 
discussed, viz., inflation 9 and strings 10 . Recently, it was pointed out 11 
that one could obtain another scenario in which inflation is followed by 
string production. 

The inflationary phase is associated with the transition from S0(10) to 
SU(3) x SU(2)|_ x SU(2)p x U(1) B _ L . It can be implemented by generalizing the 
arguments of ref. (12) where the SU(5) model is discussed. The breaking of 
B-L and the U(l) symmetry can occur during, or at the end of the inflationary 
era. The spectrum of density fluctuations produced in this scenario is 
essentially of the Harrision-Zeldovich 9 type. 

According to recent observations 13 , the value for ft obtained on scales up 
to ~ 10 15 M 0 is = 0.2 + 0.1, considerably less than unity, the value predicted 
by the new inflationary cosmology. As a reasonable upper limit for ft $c of 
superclusters 14 , we may take ft $c < 0.5. Therefore, since axions and baryons 
cluster on scales smaller than rich clusters and superclusters, their 
contribution to ft must be < 0.5. The balance of the total ft in the universe 
must therefore be in the mass density of a neutrino component which is not 
traced by the galaxy distribution if we are to have 8*1. 

We must therefore require that the neutrinos be light enough so that they 
will not cluster on scales below - 10 16 Mq. In order to arrange this, 
especially since the neutrino Jeans mass drops significantly between the 
redshift z nr when the neutrinos become nonrelati vistic and the present time, 
we invoke neutrino phase space limits using the arguments of Tremaine and 
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Gunn 15 in reverse to get an upper limit on m v . These authors find that for 
neutrinos to be able to cluster on the scale of rich clusters, their mass must 
be greater than - 4 hj^ 2 e V (where hso is the Hubble constant in units of 50 
kms"lMpc"l) . 

O p p 

The neutrino contribution to ft is ft y = 4.56 x 10 m v (eV)N^ hjlg T^ g 
where Nf is the number of neutrino flavors of approximately equal mass and 
Tg # 8 is the present temperature of the cosmic blackbody radiation in units of 
2.8 K. We require ^ to be > 0.5 so that the total 0=1. For this, one 
needs at least three flavors of neutrinos, each of approximately 3-4 eV. As 
discussed above, this situation is readily obtained in the S0(10) model. (If 
the neutrino clustering is inefficient (see discussion in Ref. 16), m v could 
be larger and Nf smaller.) 

The maximum neutrino Jeans mass for three neutrinos of roughly equal mass 
is 17 I'fj v = 2.7 x 10 18 [m v (eV) ]" 2 M 0 , which, for Nf = 3 and m y = 3.6 eV gives M*j, 
2 x 10 17 Mq. The corresponding spatial scale at present for pancaking 
structure would be ~ 150 Mpc. It is interesting to note that this scale may 
correspond to the tentative "superpancaking" scale proposed recently by 
Dekel 18 in order to attempt to account for the correlation function of 
clustering of superclusters. 19 Structure on this scale would have to 
correspond to density perturbations <5 = 6p/p just becoming nonlinear 
(^ = 0.5-1) at the present time. 

The spectrum of perturbations in a universe dominated by axions and 

neutrinos is readily estimated by adopting the arguments previously given for 

a baryon-neutrino universe 20 . It is convenient to define 5 = ft /(ft +ft ) 

a a v' 

such that 5 _< 1/2 (We assume, for simplicity, that << ft a , ft y ). 

For z < z 6 q - 0.93 x 10 4 (1-C)” 1 ft y h^ 0 T^g the neutrino Jeans mass 
decreases as (l+z) 3/2 . (Here z eq is the redshift corresponding to equal 
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matter and radiation densities in the universe.) Neutrino perturbations on 
scales below M* v are erased at z 2 z eq . The axion perturbations, however, 
grow like 


6p a 


p 

a 


6 

a 


( 1 + 2 ) 


-3 a/ 2 


(3) 


where a = (/I+ZTS -l)/6. (The growing mode solution is similar to that 
obtained for the baryon-neutrino hybrid scenerio after decoupling 20 .) Thus, 


6 

a 


(z) - * 




1+z eq,3a/2 

“T+z - ' 


(4) 


This continues until z - z^ when the neutrino Jeans mass becomes = M, 


(1+Z M J = (1+Z eq ) 


(5) 


For z < z M the overall density fluctuation 6p/p * t 2 ^ 3 « (1+z)" 1 . Thus, 


<$P . _ \ ~ e x / \ / 1+Z M\ — ext \ / 1+z eq* / M \ (2/3-a) 

- ( Z<Z M> = 5 V Z M> W 1 = 5 V z eq> ( TfT> ^ (6) 


As a rough approximation, <$ a (z eq ) " constant when M < M* y for a Zeldovich 
spectrum. (See, however, footnote 21). This gives 


— * m( 2 / 3 -“) (M < M* ) (inflation alone) (7) 

P J V 


which is an increasing function of M since a < 2/3. For M > M* y , the neutrino 
perturbations are not damped and <$p/p * M -2 ^ 3 . 

From this discussion it appears that even in the most optimistic case 
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where 5 = 1/2, a = 0.43, so that the scales between the present neutrino Jeans 
mass and M* v may not collapse before M^ v . does. We thus run into the timing 
problems which are becoming well known for the neutrino pancaking scenario. 

In particular, it is hard to envision the development of quasars 22 and 
substructure 23 with such a model, although the situation here is not as 
difficult as that with pure neutrino pancakes owing to the presence of axions 21 , as 
we discuss below. 

The presence of strings, which provide an additional source of density 
fluctuations, can eliminate the above difficulty 24 . Assume that topologically 
stable strings, with mass per unit length characterized by a superheavy (GUT) 
scale, appear at or near the end of the inflationary phase. A specific 
example showing how this could occur is shown in Ref. 11-. In the present case 
this is readily achieved either by appending a new spontaneously broken global 
U(l) symmetry to the S0(10) model or using an Eg model. Owing to the presence 
of strings, and, in particular of closed loops, 25 5 a (z e q) * M ~^ 3 for (M < M j v )* 
Substitution in eq. (6) then gives 

-jp « M (l/3-<*) (M < M^ v ) (string loops) (8) 

as compared with the results of eq. (7) when loops are not present. 

Using eq. (8) with C = 1/2, and a = 0.43 we find <5p/p * Therefore, 

if 6p/p ~ 0(1) on scales ~ 10 16 - 10 17 Mq at z=0 as suggested by Dekel 8 , scales ~ 
10 10 Mq went non-linear at z a 4, corresponding to the epoch of quasar formation. 
Thus, in the presence of axions and neutrinos, an inflationary scenario supplemented 
by strings (or wall -string systems 24 ) appears to offer a better prospect of 
explaining the observed large scale structure in the universe than one without 
strings. Of course, more detailed numerical calculations and clustering simulations 
should be performed to test this conclusion. In fact, growth of axion perturbations 
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during the radiation era ^3 wl n have the effect of increasing a to = a + e. 

This effect may be enough to make the spectrum in the case of inflation without 
strings flat at low M. In the string-inflation scenario, this effect eases the 
requirement on 8 a needed for an acceptable a e ff, making 8 < 0.5 (as indicated by the 
observations) acceptable. 

Finally, let us discuss the scenario in which we dispense with inflation and 
density fluctuations are produced solely by strings. In this case, since the 
density parameter fl need not be unity, £ can be greater than 1/2 and a can be > 

0.434. (Of course, we need have only one v flavor in the eV mass range to get 
Dekel's 18 scale.) In particular for fl a >> ft v , a = 2/3. A natural extension of 
S0(10) which gives the desired strings 25 is provided by the following breaking of Eg 
(once again the global 11(1) Peccei -Quinn symmetry is broken at the same scale as B- 
L) 

E. - S0(10) x Z, > SU(3) x SU(2). x SU(2). x U(1) R . x Z. (9) 

6 ~10 16 GeV 2 L R B-L 2 

SU(3) x SU(2 ) l U(l) x Z 2 

For Eg symmetry breaking at a scale n ~ 10^ GeV, the energy per unit length of 

2 32 2 

the strings formed is u - n » 10 GeV . (A similar result can be obtai7ed naturally 
in a Kaluza-Klein model leading to So(10) (Wetterich, private communication).) With 
this value of u, it follows from the discussion of Ref. 25 that in this scenario 
neutrino perturbations would be on the verge of becoming non-linear at the 
“superpancake" scale at the present time, as suggested by observations 18 * 19 . 

To conclude, significant axion and neutrino energy densities arise naturally in 
a class of grand unified theories. An axion-neutrmo dominated universe model for 
the formation of large scale structure may avoid the problems associated with the 
pure neutrino dominated pancake models. These models also allow for structure on 
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scales greater than that given by the pure hierarchical clustering models of galaxy 
formation, which may be desirable in view of some recent analyses suggesting the 
clustering of clusters. Finally, the prediction of the new inflationary cosmology 
that be unity can be reconciled with the observation 8 <1 in this framework, 

SC 

particularly if string loops (or string-wall systems) are present 26 . 
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